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Abstract 

Correlation mixtures of elliptical copulas arise when the correlation parameter 
is driven itself by a latent random process. For such copulas, both penultimate 
and asymptotic tail dependence are much larger than for ordinary elliptical cop- 
ulas with the same unconditional correlation. Furthermore, for Gaussian and 
Student t-copulas, tail dependence at sub-asymptotic levels is generally larger 
than in the limit, which can have serious consequences for estimation and evalu- 
ation of extreme risk. Finally, although correlation mixtures of Gaussian copulas 
inherit the property of asymptotic independence, at the same time they fall in 
the newly defined category of near asymptotic dependence. The consequences 
of these findings for modeling are assessed by means of a simulation study and 
a case study involving financial time series. 
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1. Introduction 

It is a stylized fact that financial data such as stock or exchange rate returns 
exhibit a sizeable amount of tail dependence. For that reason, student f-copulas 
with low degrees of freedom are often used to model dependence for such data. 
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Furthermore, in many instances cross-sectional correlations in multivariate fi- 
nancial time series have been observed to vary over time . Consequently, 
in some recent studies the dependence between financial variables has been 
modeled via copulas whose parameters vary themselves according to a latent 
random process The unconditional copula is then a mixture over the 
underlying parametric family according to a certain probability distribution on 
the parameter. 

In this paper we focus on the specific case of tails of correlation mixtures 
of elliptical copulas, and more specifically of Gaussian and t-copulas. This 
situation arises for instance in certain multivariate stochastic volatility models 
Q and in the Stochastic Correlation Auto- Regressive (SCAR) model of [4]. In 
the latter model, the latent cross-sectional correlation pt at time t is described 

by 

exp(27t) -I- 1 

where et, t ^'L^ are independent standard normal variables, < 1, and the 
inverse Fisher transform is chosen to keep pt in (—1, 1) at all times. This spec- 
ification is intuitively reasonable, analytically tractable and has been found to 
provide an excellent fit to financial data when used in conjunction with a Gaus- 
sian copula. Correlation mixtures over elliptical copulas being not necessarily 
elliptical anymore, the setting in our paper is not covered by the literature on 



paper is not cove: 

i fly HQ. 



tail behavior of elliptical distributions 

Let C be a bivariate copula and let (?7, V) be a random pair with distribution 
function C. Dependence in the tails can be measured by 

\(u) = u-^ C{u, u) = P(y <u\U <u), < ii 1, 

and its limit A = Imiuio X{u), the coefficient of tail dependence. (Elliptical 
copulas being symmetric, it suffices to consider lower tail dependence.) We call 
A(u) the coefficient of penultimate tail dependence and we will argue that A(u) 
may be more informative than A when their difference is large. The distinction 
between tail dependence at asymptotic and subasyrnptotic levels has already 
been made several times in the literature [ii,Q,Q,l3. 



For t-copulas with correlations generated by (|l.ip . A(m) is plotted in Figure 
[T]for values of u that go deep into the joint tail. The unconditional correlation 
is equal to E[pt] = 0.5 for all cases. The figure illustrates our two main findings 
quite well: 
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Figure 1: Coefficient of penultimate tail dependence A(n) = u ^ C{u,u) for correlation mix- 
tures of t-copulas. Mind the logarithmic scale for u. 

1. Allowing for random correlation greatly increases both X(u) and A. 

2. The speed of convergence of A(w) to A may be extremely slow, implying 
that tail dependence at penultimate levels may be significantly stronger 
than in the limit. 

Similar effects can be observed for correlation mixtures of Gaussian copulas, for 
which A(u) is much larger than its limit even at levels u that are several orders 
of magnitude smaller than the ones relevant for practice. Our findings under- 
line the importance of proper modeling of correlation dynamics and of taking 
threshold sensitivity of tail dependence measures into account. A surprising 
consequence is that Gaussian copulas with time-varying correlations do a great 
job in capturing tail dependence at subasymptotic thresholds. 

The rest of the paper is structured as follows. In Section [2] the impact of 
random variation in the correlation parameter on tail dependence measures of 
elliptical copulas is investigated. The (lack of) speed of convergence of X{u) to 
A for correlation mixtures of Gaussian and t-copulas is studied in Section [31 In 
the Gaussian case, the rate of convergence to zero is so slow that it is actually 
appropriate to speak of near asymptotic dependence. Consequences of our results 
for modeling are discussed in Section |4] through a Monte Carlo study and a case 
study. Finally, Section [5] concludes. Proofs are relegated to the Appendix. 



2. The impact of correlation dynamics 

Let {X, Y) be a (standardized) bivariate elliptical random vector, that is, 

(X,r)^ (5l,p5i + (l-p2)l/2^2) (2.1) 

where — 1 ^ p ^ 1 and where 

(S'i,S'2) = (i?cose,i?sine), (2.2) 

the random variables R and Q being independent with R > and Q uniformly 
distributed on (— 7r,7r); see for instance 0, Consider the distribution 

of the radius R as fixed. Let Cp be the copula of {X, Y) seen as parameterized 
by the correlation parameter p, that is, 

Cp{u, v) =P{U^u,V^ v), (u, v) e [0, 

with U = F{X) and V = F(Y), where F is the common marginal distribution 
function of X, Y, Si and 5*2. Let 

Xp{u) = u-^Cp{u,u) = P{U <u\V <u), 0<us^l, (2.3) 

denote the coefficient of penultimate tail dependence and let Ap denote its limit 
as u I 0, provided it exists. 

As in the SCAR model suppose now that p in (|2.ip is itself a (latent) 

random variable, independent of and 5*2. The copula C of (X,Y) is then a 
p-mixture of the copulas Cp, 

C{u,v) = j Cp{u,v)^i{dp), (2.4) 

with /i equal to the probability distribution of p on [—1, 1]. The penultimate 
coefficient of tail dependence of C is simply given by 

A(m) = C(u, u) = j u^^ Cp{u,u) p,{dp) ^ J Xp{u)p{dp), < u ^ 1. 

(2.5) 

Let the average correlation parameter be denoted by 

p = J pp.{dp) 

Provided second moments exist, p is equal to the (unconditional) correlation 
between X and Y. If for some u g (0, 1] the function p Ap(u) is convex, then 
by Jensen's inequality, 

X{u) > \p{u) (2.6) 
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whatever the mixing distribution /i. As a consequence, when correlations are 
themselves driven by a latent random process, tail dependence may be larger 
than what one may expect. Taking limits as u J, shows that this reasoning 
applies to the coefficient of tail dependence A = j\^p ^^{dp) as well. 
It remains to investigate the convexity of Xp{u) in p. 

Proposition 1. Let Xp{u) be as in (|2.3p . For every u G (0, 1/2], the function 
p I—!- Xp{u) is convex m p G [0, 1]. 

Graphs of Xp{u) for fixed u (not shown) suggest that the restriction to p & 
[0, 1] cannot be avoided. Still, under some mild conditions on the distribution 
of R, the convexity is actually true on p e [—1, 1] for all sufficiently small u > 0. 

Proposition 2. Suppose that the distribution function Fj^ of R in (j2.2p has 
unbounded support and is absolutely continuous with density fn. If 

liminf^i^>l, (2.7) 

r^oo Fji{r) 

then there exists Uq G (0,1/2) such that the function p i— s- Xp{u) is convex in 
p G [—1, 1] for every u G (0,uo). // the limit Ap ~ lim„j^o Ap(u) exists, then Ap 
is convex m p G [—1, 1]. 



Condition (|2.7p is verified as soon as the radial density function is regu- 
larly varying at infinity of index —a — 1 for some a > 1, in which case the liminf 
is actually a limit and is equal to a. For the bivariate t-distribution with > 
degrees of freedom, the radial density is given by /i?(r) — r (1 + r^/i^)^"/^^^ 
for r > 0. Since fnir) = (c + o{l))r~'^~^ as r ^ oo for some constant c > 0, 
we find r ffj{r) / F ii{r) ^ as r ^ oo. Hence, condition (|2.7p is satisfied if 
V > 1. For the bivariate Gaussian distribution, the radial density is given by 
fR{^) — ^ exp(— r^/2) for r > 0, which implies r f ji[r) / F f>{r) — ^ oo as r — > cxi, 
so that (|2.7p is satisfied again. 

The difference between \{u) and Ap(u) can be seen in Figure [H whereas 
Figure [2] illustrates it for the coefficient of tail dependence A. For the (static) 
i-copula with v > degrees of freedom and correlation parameter p G (—1, 1), 



the coefficient of tail dependence is given in [20|, |2l|, [2^ by 



A.^p = 2t,+i(-x^7TTv/T^/Vl + p), (2-8) 

where t^ denotes the cumulative distribution function of the <-distribution with 
p > degrees of freedom. In Figure [21 Student f-copulas are considered with 
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Figure 2: Coefficients of tail dependence of static and dynamic t-copulas. 

stochastic and static correlations such that both have an unconditional corre- 
lation coefRcient oi jj ~ 0.5. The degrees-of-freedom parameter i' varies from 2 
to 30. Two different correlation dynamics are considered: the SCAR process in 
(jl.ip . and the uniform distribution on (0, 1). The dynamic models clearly lead 
to significantly higher tail dependence than the static ones. 

3. Penultimate tail dependence 

So far we have shown that correlation mixtures over elliptical copulas exhibit 
stronger tail dependence than static ones with the same unconditional correla- 
tion. For practical purposes, however, the penultimate coefficient X{u) is more 
suited to assess the risks of joint extremes than its limit A, sec for instance 
In this section we study the rate of convergence of A(w) to A for corre- 
lation mixtures of the Student i-copula (Section 13. ip and the Gaussian copula 
(Section [321) ■ 

3.1. Correlation mixtures of t-copulas 

Let v > and —1 < p < 1. Let Zi, Z2, S be independent random variables 
such that the Zi are standard normal and vS^ has a chi-square distribution with 
v degrees of freedom. Put 

X = S-^Z,, Y = S-^pZ, + ^/T^Z2). 
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The distribution of the random vector {X, Y) is bivariate t with ly degrees of 
freedom and correlation parameter p. Its copula, C, is the bivariate t-copula 
with parameters u and p. 

Now suppose as in Section [2] that p is itself a random variable with range in 
(—1,1) and independent of Zi, Z2, S. The unconditional copula, C, of {X, Y) is 
then a correlation mixture of bivariate f-copulas with fixed degrees-of-freedom 
parameter ly. 

For nonrandom p £ (—1,1), the coefficient of tail dependence X^ p is given in 
(|2.8p . For general, random p, we have A — E[Ai/,p], the expectation being with 
respect to p. 

Figure [1] suggests that the rate of convergence of X{u) = u^^ C{u,u) to A 
may be slow, especially when v is large. This is confirmed in the following result. 
For real x, let — max(a:;, 0) be its positive part. 

Proposition 3. Let C be a correlation mixture of bivariate t-copulas with 
degrees-of-freedom parameter 1/ > 0. We have 

X{u) = u-^ C(u, u) = A + 7 u^/" + 0(^2/"), u i 0, 

where j is a positive constant depending on v and the distribution of p given in 
(|A.3P below. 

The rate of convergence of A(m) to its limit A is of the order 0{u^^'^). The 
larger the slower this rate. Since 7 is positive, X{u) may therefore be (much) 
larger than A. As a result, at finite thresholds, the tail may look much heavier 
than it is in the limit. This in turn may cause estimators of v to be negatively 
biased. 

3.2. Correlation mixtures of Gaussian copulas 

Let X, Z^p be independent random variables, X and Z being standard nor- 
mal and p taking values in (—1, 1). Put 

Y = pX + p^f/^Z. (3.1) 

Conditionally on p, the distribution of {X, Y) is bivariate normal with zero 
means, unit variances, and correlation p. The unconditional correlation of X 
and Y is corr( JT, F) = E[p] = p, and the copula of {X,Y) is a correlation 
mixture of Gaussian copulas: 

C{u,v) = F.[Cp[u,v)], (u,i;) e [0,1]2, (3.2) 
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the expectation being with respect to p, and with Cp denoting the bivariate 
Gaussian copula with correlation p. As before, put A(u) — u^^ C{u,u) and 
Xp{u) = Cp{u, u). 

By the assumption that p < 1 almost surely and since the coefficient of 
tail dependence of Cp is equal to = lim„|o Ap(m) — 0, the coefficient of tail 
dependence of C is A = E[Ap] too. So just like Gaussian copulas, correlation 
mixtures of Gaussian copulas have asymptotically independent tails. 

To measure the deg ree of tail association in case of asymptotic independence, 
Ledford and Tawn fl6l| introduced the coefficient 

Mio logC(w, u) 

the existence of the limit being an assumption. In this case, one can write 

C(u,w) = ui/"L(u), (3.4) 

the function L satisfying lim„|o log L(u) / log u = 0, that is, for all e > we have 

lim u'^ L(u) = 0, limu~^ L(u) = oo. 

uio uio 

Typically, the function L is slowly varying at zero: lim„, m L(ux) / L(u) = 1 for all 
a; > 0. The coefficient r] is related to the coefficient x in UM through x = 2?7 — 1. 

The pair of coefficients (A, rj) measures the amount of tail dependence across 
the classes of asymptotic dependence and independence. The following two 
cases are most common: 

1. Asymptotic dependence: If A > 0, then necessarily rj — 1. 

2. Asymptotic independence: If 77 < 1, then A = 0. 

For the bivariate Gaussian copula with correlation p < 1, for instance, we 
have A = and rj = {p + l)/2, so x — P- One may then be tempted to believe 
that for correlation mixtures of Gaussian copulas the coefficient x is equal to 
p. However, this is false. Instead, a new situation is encountered, one that is in 
between the two cases described in the preceding paragraph. 



Proposition 4. Let C he as in (j3.2p . If p < 1 almost surely but the upper 
endpoint of the distribution of p is equal to 1, then for every £ > 0, 

limu^"^ X(u) = 00. 

As a consequence, although we have A = (asymptotic independence), the 
Ledford-Tawn index is equal to rj ^ 1; likewise, x = 1- 
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Figure 3: The coefficients X{u) and r] for correlation mixtures of the Gaussian copula. 

According to Proposition [31 the rate of decay of \{u) to A = is slower 
than any positive power of u. So even though the tails of C are asymptotically 
independent, we are as close as one can get to the case of tail dependence. 
Therefore, we coin this type of tail behavior near asymptotic dependence. As 
far as we know, this situation has not yet been encountered in the literature; 



see for instance the extensive list of examples in 15|. Note that the function 
L{u) in (|3.4[) is equal to \{u). 

By way of example, the left-hand panel in Figure [3] shows \{u) for tail 
probabilities u ranging from 10"'^ to 10^^", the correlation parameter being 
driven by the SCAR model (jl.ip for various choices of the parameters, but 
always with p — 0.5. Especially when the variation in p is large, the coefficient 
\{u) remains quite sizeable even very far in the tail. To investigate the speed 
of decline of A(u) numerically, we write 

log \{u) = (3 log u + log L[u) 

with [3 = l/rj — 1 and we treat log L (it) as a constant, its rate of change being 
much slower than the one of logw. The slope (3 can be estimated by a simple 
least squares regression, from which 77 can be recovered. We estimated 77 in this 
way by making u vary over intervals [ul, uu] — [u^^ , u^^"^] for k ranging from 
3 to 10 with steps of size 0.01. The corresponding estimates of 77 as a function of 
ul are depicted in the right-hand panel Figure [31 For the constant correlation 
case we find indeed 77 ~ (1 -t- /9)/2 = 0.75. In contrast, when correlation is not 
constant, 77 seems to be converging to 1 as it | 0. 
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4. Consequences for modeling 

In this section we illustrate the consequences of our results for modeling 
data with static versus dynamic t-copulas. Of particular interest is the effect on 
inference on tail dependence when correlation is falsely assumed to be constant. 
To fit a static t-copula, we estimate v by maximum likelihood and p by the 
method of moments, the likelihood function for v being given by 

L{v I u,v)^U{t-^{u), tz\v),p), (4.1) 

where t denotes the density of the bivariate ^-distribution and p is the empirical 
correlation of the pairs (t^^(J7i), t^^(Vi)). 

4-1. Simulation study 

We simulated data from a i-copula with a fixed degrees-of- freedom parameter 
V but stochastic correlations. The correlations were drawn from the SCAR 
process in (jl.ip . We chose (3 ~ 0.97 and a — 0.05,0.1,0.15,0.2, whereas a 
was chosen such that the unconditional correlation was in all cases equal to 
p — 0.5. The values for the degrees-of-freedom parameter v were 5, 10, 50, and 
cx), the last case yielding the Gaussian copula. The sample size and the number 
of Monte Carlo replications were chosen to be equal to 1000. For each of the 
1000 simulated data sets we fitted a static t-copula and computed both the 
implied coefficient of tail dependence Xo^p as well as the one of penultimate tail 
dependence Ap_p(u) at u = 0.01, a tail probability of potential practical interest. 

The bias for the estimates of A and \{u) is reported in Table [TJ Overall, 
X{u) is estimated much more accurately than A and both are estimated best 
when V is low. The bias generally increases in v and for A it also increases 
in the volatility a of the correlation process. Notably, the estimates suggest 
quite strong tail dependence even when the true conditional copula is Gaussian. 
The estimates of p and v (not shown) show that the positive bias in A is due 
a negative bias in whereas the overall correlation p — 0.5 is only slightly 
underestimated with average estimates ranging from 0.48 to 0.5. It seems as 
though the estimate of the degrees-of-freedom parameter v is such that the 
implied coefficient of penultimate tail dependence A(O.Ol) matches the true one 
and that the consequence of this is a severe underestimation of the coefficient 
of tail dependence A. 

These findings continued to hold for other sample sizes and other correlation- 
driving processes. 
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Table 1: Bias estimating A(u) and A using t-copulas 



A(n) 


a = 0.05 


a = 0.1 


a = 0.15 


a = 0.2 


v= 5 


-0.011 


-0.001 


0.002 


-0.022 


v = lQ 


-0.027 


0.014 


0.022 


0.008 


v = 20 


-0.082 


0.005 


0.029 


0.031 


l> = oo 


-0.096 


0.000 


0.034 


0.030 


X 










v= 5 


0.002 


0.029 


0.029 


0.012 


J/ = 10 


0.027 


0.081 


0.107 


0.097 


J/ = 20 


0.039 


0.142 


0.213 


0.235 


v = oo 


0.022 


0.131 


0.240 


0.303 



Note. This table reports on the bias for A and A(O.Ol) for i-copulas 
when the correlation is assumed to be constant when they are in fact 
generated by the SCAR model 11.11 1. The sample size and the number 
of Monte Carlo replications are equal to 1000. 



^.2. Application to stock market and exchange rate returns 

We consider international stock market returns at daily and monthly fre- 
quencies and daily exchange rate returns. The data sets are daily returns of the 
Dow Jones industrial index (DJ) and the NASDAQ composite index (NQ) from 
March 26, 1990 until March 23, 2000, daily returns of the MSCI index for France 
(Fra) and Germany (Ger) from October 6, 1989 until October 17, 2008, monthly 
returns of the Datastream stock index for Germany, Japan (Jap), the UK and 
the US from January 1974 until May 2008, and daily exchange rate returns of 
the Euro (EUR), British pound (GBP) and Japanes Yen (JPY) against the US 
dollar from January 1, 2005 until December 31, 2008. Returns are calculated as 
100 times the first difference of the natural logarithm of prices. 

To model the marginal distributions, we opt for the stochastic volatility 
model 23|, |2j| because of its natural connection to the SCAR model for the 
correlation dynamics (jl.ip . The basic stochastic volatility model for the return 
rt at time t = 1, . . . , T is given by 

n = exp(/it/2) et, ht = ah + h h-i +Tr]t, 

where et and r]t are independent standard normal random variables, uncorrelated 
with the innovations driving the dependence process (jl.ip . Estimation of the 
model is done by simulated maximum likelihood using the efficient importance 
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sampler 25|. The static i-copulas are fitted as described above, whereas for 
the time-varying model we condition the estimate of v on the correlations pt 
estimated using the SCAR model, yielding estimates pt and v\pt respectively. 

We also compute the implied coefficient of tail dependence A and the coef- 
ficient of penultimate tail dependence A(u) at various economically interesting 
levels. In particular, the levels we consider correspond to exceedances that are 
expected to occur once a year, once a decade and once a century. We denote 
the corresponding estimates by Ayearj Adec and Acent- These can be interpreted 
as the probabilities that one market makes a certain large downward movement 
conditionally on the other market doing the same. The coefficient of tail de- 
pendence A, on the other hand, denotes the probability of one market crashing 
completely, i.e. dropping to a level of zero, conditional on the same event for 
the other market. The latter scenario is economically rather unrealistic and is, 
in our view, relevant only for individual stocks, not for complete markets. 

The results can be found at the end in Table IA.2I The most striking find- 
ing is that the estimated degrees-of-freedom parameter v is significantly larger 
when correlations are allowed to be stochastic, and in many cases even virtu- 
ally infinity, corresponding to the Gaussian copula (we report oo whenever the 
upper bound of 400 in the optimization routine was obtained). This implies 
that a large part of the fat-tailedness can be captured by random correlations. 
For the estimated measures of tail dependence, two things are notable. First, 
the penultimate version A(m) is much larger than the limiting one A. This sug- 
gests that A may be a too optimistic measure for assessing the risk of spillovers 
of large downward movements across financial markets. Second, although the 
fitted models that allow for time- varying correlations have a lower limiting coef- 
ficient of tail dependence than the static ones, at practically relevant quantiles 
these models capture the dependence in the tails of the distribution quite well. 
So just like in the simulation above, the static and the dynamic models both do 
a good job in matching penultimate tail dependence, but the static t-copulas 
may lead to overestimation of the coefficient of tail dependence A. 



5. Conclusions 

We have studied tail dependence properties of correlation mixtures of el- 
liptical copulas, a situation which occurs when the correlation parameter is 
itself driven by a latent random process. We have shown that the coefficient 
of (penultimate) tail dependence is larger than for ordinary elliptical copulas 
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with the same unconditional correlation. Furthermore, for Gaussian copulas 
and i-copulas, tail dependence at sub-asymptotic levels quantiles can be sub- 
stantially larger than in the limit. In a simulation study we found that ignoring 
the dynamic nature of correlations when estimating i-copulas leads to biased 
estimates of the coefHcient of tail dependence. Our empirical application showed 
that estimates of the degrees-of-freedom parameter of a i-copula are much lower 
when assuming a static correlation than when conditioning on dynamic corre- 
lations. At the same time the models based on dynamic correlations produce 
similar dependence in the tails at economically relevant quantiles, but lower tail 
dependence in the limit. 

A notable discovery was that under some fairly weak conditions, correlation 
mixtures of Gaussian copulas have tails that are so close to being asymptotically 
dependent that they give rise to the newly defined category of near asymptotic 
dependence. In practice this implies that it is virtually impossible to distinguish 
such copulas from ones that have asymptotically dependent tails such as t- 
copulas. 

These findings suggest that for practical purposes the Gaussian copula is 
more attractive than as often stated in the literature, as long as one accoimts 
for the (empirically observed) fact that correlations vary over time. This finding 
can be seen as an analogue to the effect that conditionally Gaussian models for 
time-varying volatility such as GARCH and stochastic volatility models can 
create heavy tails in the margins. Thus, conditionally Gaussian models are 
more than just a simplifying approximation in a multivariate setting, as they 
are able to capture the tails both in the margins as in the copula. 
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Appendix A. Proofs 

Proof of Proposition]^ Write p = sin 7 G [^1,1] for 7 = arcsinp e 
Then (1 - p^jVa = C0S7 e [0, 1], so that 

Y ^ pSi + {1- p^Y/^S2 = R sin 7 cos e + i? cos 7 sine = R sin(e + 7), 

yielding the representation 

(X, Y) = (i? cos e, R sin(e + 7)) . 

Let i ^ be such that 1 — F{t) — F{—t) = u, where F is the marginal distribu- 
tion function of X and Y . Since the distribution of (X, Y) is symmetric around 
zero and since R and O are independent, 



/oo 
p 



cos Q > - , sin(6 + 7) > - 
r r 



dFR{r), 



with Fj^ the distribution function of R. As a consequence, it is sufficient to 
show that for fixed z S [0, 1), the function 

p P[cos0 > z, sin(6 + 7) > z] 

is convex in p G [0, 1], where 7 = arcsinp G [0, ^] and with O uniformly dis- 
tributed on (— tTjTt). Write z — cos a G [0, 1), so a = arccosz G (0,7r]. Then 

cos 9 > cos a -i^^^ —a < <d < a 

whereas 

sin(6 -f 7) > cos a = sin(^ — a) <;=^ ^ — a<Q + j< ^ + a. 
Joining these two double inequalities and using the fact that ^ ~ j ^ yields 



P[cose > z,sin(e-|-7) > z] ^ P(-Q;-f f-7 < 9 < a) = max(2Q;- ^ -f 7, 0). 



Since 7 = arcsinp is convex in p G [0, 1], the result follows. □ 

Proof of Proposition [H We keep the same notations as in the proof of Proposi- 
tion [T] Now let < u < 1/2 so that t = F-'^{1 - u) > 0. We have 



Cp(u,u) 



^ ^ max(2a - f + 7, 0) dFn{r) 
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where a = arccos(t/r) and 7 — jip) — arcsinp. Observe that 
2a-f+7>0-^r> — ^ — — = r{p) 



cos 



((f-7)/2) 



We find 



dp ^ ^ ' dp ' 

2n -^C.iu, u) = -fn{rip)) ^ ^ + Fn{r{p)) 

We have 

^ (1 _ p2-j-l/2 <Pl{p) ^ P djjp) 

dp ' dp'^ 1 — p^ dp 

Some goniomctric juggling yields cos{(^ — 7)/2} = {(1 + p)/2}^/^, whence 

dp 2 1 + p 

Writing i'{r) — r fji{r)/ F{r), we conclude that 

For p < 0, we have r{p) > t. If u is small enough so that t is large enough so that 
i'{r) > 1 for all r ^ t, then the factor between big brackets on the right-hand 
side of the last display is positive for all p € (—1, 1). Hence Cp{u,u) is convex 
in p. □ 

Proof of Propositionl^^ Let a; > 0. We have 



P{X >x,Y >x)= F{S-^Zi > X, S'^ipZi + v/1 - p^Z2) > x) 
= P(min(Zi,pZi + ^/T~^Z-2)S-^ > x) . 



Put W = min(Zi , pZi + \/\~~fp- Z2) and W+ = max( W, 0) . Then we can rewrite 
the above equation as V{X > x,Y > x^ — P(W+5'~^ > x). We need precise 
information on the upper tail of S^^: see Lemma [S] We can then proceed as 
follows: 

V{X >x,Y >x)= [ P{S-^ > x/w)dP{W < w) 
{x/w)-" (1 - 6„ {x/w)-^) dP{W w) 

(0,00) 

A(x/w)dP(Vl^ < w) 

(0,00) 

Wl]~Kx-'' E[Ty^+2]) + I 

(0,00) 



x""" (E[M^^] - 6^ x-"^ E[Ty|+2]) + j A(a;/w) dP(Vl^ w). 

J(O.oo) 
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The remainder term is 



0^ / A{x/w)dP{W s^w) s^c^x^"-^ E[W:;.+^]. 

J(0,oo) 

We find that, as x oo, 

P{X > x,Y > x) =a^ x-" (E[W]:] - x~^ E[Ty|+2] + ©(x"*)) . 

The marginal tail of X can be represented in the same way: it suffices to 
replace W in the preceding display by a standard normal random variable Z , 
that is, as a; — !■ cx), 

P(X > a;) = V{S-^Z > x) = x^" (E[Z^] - K E[Z;+2] + 0{x~^)). 

It follows that as a; — !■ cx), 

_ Y.[Wl] ~ K x-"^ E[PF;+^] + 0{x-^) 
~ E[Zl]-b,x~^ E[Zl+^] + 0{x-^) 

- E[Z^] +'^" (E[z^])^ ^""-'^ 

As a consequence, the coefficient of tail dependence is given by 

E[A.,p] = A = lim PiY>x\X>x)^ -±-±(. (A.2) 

X^OO Cj[Z^\ 

Moreover, the expansion in (|A.ip gives us a handle on the rate of convergence 
of A(m) towards A. Let be the quantile function of the i-distribution with v 
degrees of freedom. By symmetry of the upper and lower tails, 

A(?i) = p(y > tz^{i - m) I a: > t-^{i - u)) 

_^ , ^ E[iy-]E[Z-+^]-E[W^;+^]E[Z-] 2 

+ 0({t~\l-u))-') 
Since 1 — (x) = P(A" > a;) x^'^ E[Z^] as a; — > oo, we have 
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(By /(y) ~ g{y) we mean that f{y)/g{y) ^1-) We obtain \{u) = A + 7 u^/'' + 
oiu^^") with 

7 = 0^ ; 72 1^1/ ) 

(E[Z^])' 

Equating fO]) and jM]) (for non-random p) we find 

E[VK;] =E[Z;] E[A,,p], 

the latter expectation being with respect to the random variable p. As a conse- 
quence, 

E[W-] E[Z;+^] - E[W;+^] E[Z-] _ E[Z;+^] 

the sign of which is positive, for \y p is decreasing in v. Furthermore, by (jA.4[) 
below and the identity Y[z + 1) = zY{z\ 

E[zr^] ^r((. + 3)/2) ^, + 1 

E[z;] -z^r((. + i)/2) 2 ^+ ■ 

We obtain 

7= (a, E[z;])"'/"6,(i. + l) E[A,,p-A,+2,p]. (A.3) 
with ay and given in (|A.5p and 

E[Z;] = -^y^ ,^e-^/2dz=-^r((z.-M)/2). (A.4) 

□ 

Lemma 5. for > 0, Zet S be a positive random variable such that vS"^ has a 
chi-square distribution with degrees of freedom v. Then for all z > 0, 

P(5-i > z) = ay - by z-"^) + A(2) 

where 



r(i./2 + l)' "''"i./2 + l 

and A(z) s$ z"'^--* /or some positive constant 
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Proof. Since vS'^ is chi-squared with u degrees of freedom. 

— 2 

1 



"'0 



r(i//2)2'^/2 



For z ^ 0, we have e ^ — {I — z) = {1 — e y)dy = e ^ dxdy and thus 

e-^ - (1 - z) ^ z72. We find that 

^ \ Z"^' y-/2-i = constant z^"-^ 



the constant depending on v. The integral can be computed as follows: 



1 



r(i./2)2-/2 7o 



y^'^-'{l-y/2)dy 

1 ({vz-^y/^ l{vz-^Y/^+^' 



r(j//2) 2''/2 V v/2 2 v/2 + l 

with Oi, and hi, as in (|A.5p . □ 

Proof of Proposition^ Fix e > 0. We have to prove that u^^^ = o[C{u,u)) as 
u 10. RecaU the representation {X,Y) ^ {X,pX + p^f/^Z) in Fix 
< po < 1- By symmetry, if < u < 1/2, 



C(u,u) = P[X > $"1(1 -u), pX+^Jl- p^Z > $"1(1 - u)] 

^ P[X > $-1(1 - u), po^ > *"^(1 -u), p> Po, Z>0] 

= ip(p>Po) P[X>$-i(l-u)/po]. 

Let ip{x) = (27r)~i/2 exp(— a;2/2) denote the standard normal density function 
and write $ = 1 — $. From Mill's ratio, $(a;) = (l + o(l)) x~i <p(x) as x — > oo, 
it follows that $"1(1 - u) = (l + o(l)) (-21ogw)i/2 as u I and therefore 

P[X > $-1(1 - u)/po] = $($-1(1 - u)/po) 

^(i+o(i))^(^";/^-"y^°) 

V ^ $-i(1-m)/po 
= (1 + o(l))(2^)-i/2po(-21ogu)~i/2u(i+°(i))/''o. 

If Po G (0, 1) is chosen such that 1/pq < 1 + e, then P[X > $-^(1 - u)/po] is 
indeed of larger order than u^^^ , as required. □ 
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Table A. 2: Finite level and asymptotic tail dependence of static and dynamic t-copulas 



Static correlations Dynamic correlations 



Data 


V 


'^year 




'^cent 


A 




Ayear 


Adec 


•^cent 


A 


Daily stock market 






















DJ-NQ 


6.06 


0.36 


0.33 


0.32 


0.31 


oo 


0.26 


0.18 


0.13 


0.00 


Fra-Ger 


2.89 


0.51 


0.50 


0.50 


0.50 


17.51 


0.30 


0.25 


0.22 


0.21 


Monthly stock market 






















Gcr-Jap 


12.39 


0.24 


0.11 


0.05 


0.02 


oo 


0.23 


0.08 


0.04 


0.00 


Gcr-UK 


6.22 


0.36 


0.25 


0.19 


0.18 


15.95 


0.37 


0.37 


0.36 


0.08 


Ger-US 


3.93 


0.39 


0.31 


0.28 


0.27 


11.20 


0.35 


0.23 


0.15 


0.11 


Jap-UK 


10.32 


0.26 


0.13 


0.07 


0.04 


41.03 


0.25 


0.10 


0.03 


0.00 


Jap-US 


9.91 


0.26 


0.13 


0.07 


0.05 


oo 


0.26 


0.11 


0.06 


0.00 


UK-US 


14.85 


0.41 


0.25 


0.14 


0.07 


oo 


0.47 


0.43 


0.42 


0.00 


Exchange rates 






















EUR-GBP 


8.03 


0.37 


0.33 


0.31 


0.28 


39.62 


0.31 


0.23 


0.18 


0.03 


EUR-JPY 


5.51 


0.20 


0.18 


0.17 


0.16 


oo 


0.14 


0.09 


0.06 


0.00 


GBP-JPY 


4.55 


0.17 


0.16 


0.15 


0.15 


oo 


0.21 


0.16 


0.13 


0.00 



Note. This table reports on the estimates of the degrees-of-freedom parameter v oi a. t-copula conditional on constant correlation 

(column 2) and on time- varying correlations driven by a SCAR process (column 7). The remaining columns show measures of tail 
dependence at finite and asymptotic quantiles. The margins are fitted via Gaussian stochastic volatility models. 



